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Abstract 

Bayes  estimators  are  derived  by  means  of  the  Dirichlet  process 
hyperprior  approach  for  general  empirical  Bayes  problems.  For  any 
sample  size,  these  estimators  are  expressed  concisely  as  ratios  of 
two  multidimensional  integrals.  A  numerical  example  on  Poisson 
sampling  is  given. 
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1.  Introduction 


In  the  general  setting  of  empirical  Bayes  problems,  it  is  assumed  that 
the  unobservable  parameters  { 6 ^ }  ,  i  =  1,  ....  n  ,  are  taken  indepen¬ 
dently  from  an  unknown  distribution  6  ,  and  that  associated  with  each 
6.  ,  a  random  variaole  -x.  is  observed  with  known  probability  density 
f(x-|6-)  with  respect  to  some  c-finite  measure  v  on  the  real  line.  It 
is  also  assumed,  given  the  {P^}  ,  the  observations  x^  ,  i  =  1,  ...,  n  , 
are  independent.  We  intend  to  make  inferences  about  {P.}  or  6  from 
the  observations. 

Several  approaches  to  estimating  G  or  {fc^}  are  available.  One 
approach  is  to  use  the  observed  data  to  estimate  the  mixing  distribu¬ 
tion  G  and  use  this  estimated  G  as  a  Bayes  prior.  To  estimate  this 
prior,  most  authors  assume  a  parametric  representation  of  the  prior  with 
unknown  parameters  estimated  by  the  data.  There  is  another  approach  to 
empirical  Bayes  problems,  namely  the  Dirichlet  process  hyperprior  approach, 
where  the  {e^}  are  taken  independently  from  a  random  distribution  G 
whicn  is  chosen  from  tne  Ferguson's  Dirichlet  process  (1973)  indexed  by 
a  finite  measure  a  .  The  measure  a  usually  represents  the  statisti¬ 
cian's  prior  belief  about  G  .  The  Bayes  estimator  of  G  or  {e^} 
can  be  derived.  As  pointed  out  by  Anderson  and  Louis  (1979),  this 
approach  is  potentially  superior,  since  the  construction  of  estimators 
does  not  depend  on  a  specific  form  for  the  prior.  Another  desirable 
feature  is  pointed  out  by  Berger  (1980b,  p.  83):  the  statistician  can 
combine  subjective  information  q  and  past  data  to  estimate  G  and 
(6 ^ }  ,  unlike  the  usual  empirical  Bayes  approach,  where  the  unknown 
prior  parameters  are  completely  estimated  by  the  data. 
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This  Dirichlet  hyperprior  approach  was  first  proposed  by  Antoniak 
(1974)  and  subsequently  studied  by  Berry  and  Christensen  (1979), 
Anderson  and  Louis  (1979).  The  usefulness  of  this  approach  had  been 
limited  in  the  past  bypthe  following  two  deficiencies.  1.  No  concise 
expressions  for  the  proper  Bayes  estimators  had  been  given  for  n  >  3 
due  to  the  complex  bookkeeping  and  labor  involved  in  deriving  them. 

2.  No  satisfactory  numerical  methods  had  been  developed  in  evaluating 
those  estimators.  Kith  the  work  of  Lo  (1978),  Bayes  estimators  can  be 
derived  for  any  n  .  The  purpose  of  this  note  is  to  exhibit  Bayes 
estimators  of  {6^}  for  arbitrary  n  .  It  can  be  seen  from  equation 
(1)  of  Section  2  that  each  of  the  Bayes  estimators  can  be  written  as 
a  ratio  of  n-dimensional  integrals.  These  integrals  are  hard  to 
evaluate  explicitly  due  to  the  high  dimensionality  and  the  fact  that 
the  integrands  are  peaked  in  a  small  region  of  the  parameter  space.  In 
a  recent  article,  Kuo  (1985)  proposes  to  circumvent  this  problem  by 
1.  decomposing  each  of  the  multidimensional  integrals  into  a  weighted 
average  of  products  of  one-dimensional  integrals  and  2.  approximating 
each  of  the  weighted  averages  by  an  importance  sampling  Monte  Carlo 
method.  It  is  easy  to  implement  this  method.  Moreover,  the  Monte 
Carlo  estimator  has  been  demonstrated  to  work  well  in  terms  of  effi¬ 
ciency  and  precision.  For  the  detailed  method,  statistical  analysis 
and  numerical  examples,  see  Kuo  (1935).  A  numerical  example  on  Poisson 
sampling  using  the  method  of  Kuo  is  given  here. 

See  Robbins  (1955)  for  the  pioneer  development  of  empirical  Bayes 
methods  of  estimation.  See  Susarla  (1982)  for  an  expository  article 
on  empirical  Bayes  theory  which  also  includes  some  of  the  recent 
developments  in  this  area. 
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Derivation  of  Bayes  Estimators 


To  derive  the  Bayes  rule  of  {8.}  ,  let  us  make  the  following 
assumptions: 

i.  Let  a  be  a  finite  measure  with  finite  second  moment  on  a 

measurable  space  (1R,  B)  with  1R  the  real  line  and  B  the  o-field 

•> 

of  Borel  sets.  An  unknown  distribution  G  is  chosen  from  a  Dirichlet 
process  with  parameter  a  . 

ii.  Given  G  ,  tne  unobservable  6-j,  ....  8n  are  chosen  indepen¬ 
dently  from  G  . 

iii.  Given  G  and  6  =  (8^,  ....  8n)  ,  the  observations 

n 

x  =  (x,,  x  )  have  density  f(x  |  8)  =  H  f.(x,  |  e.)  idependent 
i  n  ~  -  i  =  l  1  1  1 

of  G  ,  where  for  all  i  =  1,  ....  n  ,  f^  is  a  density  dominated  by 
v  and  f.(x  |  6)  is  measurable  in  6  for  all  x  . 

iv.  The  loss  function  is  given  by 


L(8,  a)  * 


,h  (e i  -  a‘> 


REMARK  1.  Note  that  in  Assumption  iii,  we  allow  different  f^  in  the 
model.  This  has  the  advantage  of  incorporating  individual  characteris¬ 
tics  of  the  x^'s  in  the  model,  such  as  combining  normal  and  gamma 
components  (see  Berger,  1980a),  or  treating  x^  with  different 
variances,  etc. 


REMARK  2.  Note  that  given  G  ,  the  x^  are  independently  distributed 
according  to  ff ■ (x . j 8 )G(dr)  .  This  is  essentially  the  random  density 
considered  by  Lo.  The  main  difference  is  that  his  objective  is  to 
estimate  the  mixing  distribution  G  and  various  functions  of  G  ;  our 
objective  is  to  estimate  the  e.’s  . 
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REMARK  3.  In  addition  to  lo's  work,  nonparametric  density  estimation 
has  also  been  studied  by  ferguson  (1983).  The  Monte  Carlo  method 
described  in  Section  1  was  adapted  by  Ferguson  to  compute  the  density 
estimator.  The  feasibility  of  this  method  and  error  reduction  tech¬ 
niques  were  further  illustrated. 

REMARK  4.  Empirical  Bayes  estimation  (as  opposed  to  Bayes  estimation) 
of  the  density  function  f(x[6)  =  /f(x[6)G(de)  described  in  Remark  2 
has  also  been  stuaied  by  Ghorai  and  Susarla  (1982). 


Before  proving  the  main  result,  we  first  define  some  notation  and 
state  a  lemma.  Let  a  and  G  be  defined  as  in  Assumption  i.  Then 
denotes  the  probability  measure  on  (3, A)  yielding  the  random 
distribution  G  ,  where  ®  is  the  space  of  distribution  functions  on 
(1R,8),  and  A  is  the  o-field  of  Borel  sets  in  the  Levy  metric.  The 
following  lemma  is  from  Lo  (1978,  1984). 


LEMMA  1 .  Let  6^  be  chosen  from  G  .  If  g(6-|,G)  is  a  quasi -integrable 
function  with  respect  to  the  joint  probability  G ( d9 ^ ( dG )  defined  on 
(1R  x  3,  BxA) ,  then 

a ( de , } 

(g  /m9(e,,G)G(de,)Pa(dG)  •  Jr  /B9(e,.G)Pttt£ei WG)-^  . 

We  are  now  ready  to  exhibit  the  Bayes  estimator  e  . 
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THEOREM  1 .  Given  the  Assumptions  i  through  iv,  the  Bayes  estimator  is 

A  A 

given  by  6  =  (e^,  ....  6^)  ,  where 


for  all  k  =  1 ,  . . . ,  n  . 


Proof:  Let  h(x,5)  denote  the  joint  density  of  x  and  e,  h(x) 
denote  the  marginal  density  of  x  ,  and  p^(6)  denote  the  joint 
distribution  of  6's  given  G  .  Let  P  ,  ®  be  defined  as  in  Lemma  1. 
Then  the  Bayes  estimator  under  squared  error  loss  for  6k  is  given  by 
eR  =  f(&k  |  x)  ,  where 


E(ek|x) 


!  ek  h(x,e)de 
h(x) 


;S  ;iRn  6k  f(xie)o6(de)Pa(dG) 

a”;iRn 

's  (jk  4  •  '«fkV\!ekW®kiV«i 

~s  ,5,  \  fi(x,le,:>c(df,)rQ(dG) 

/ 

/ 


by  repeated  use  of  Lemma  1. 
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REMARK  5.  It  was  shown  by  Antoniak  (1974),  the  posterior  distribution 
of  G  given  x  is  a  mixture  of  Dirichlet  processes: 

S  Ix'/Pfa*  j,  66j)dF?|x 

where  Fc ,  is  the  posterior  distribution  of  6  given  x  .  It  can  be 
e  I  x  ,  ~ 

seen  from  Theorem  1  (or  Remark  2  of  this  paper  and  Theorem  1  of  Lo  (1978)) 
that  the  posterior  mixing  distribution  of  6  may  be  written  as 


Fe;x<cl  ■ 


V-'i5,  f<(x’'8’)  J 


i-1 

&  +  I  6e 

0  =  1  eo 


(dei) 


n  ^ (xi |e.)  n 

Rn  i=l  1  1  1  i  =  l 


n  r  T-l 

a  +  I 


0  =  1 


«*i) 


for  all  C  e  En  ,  where  (lRn,En)  is  the  n-fold  product  measure  space 
of  ( 1R  ,6)  . 


3.  Numerical  Example 

To  illustrate  the  use  of  the  nonparametric  Bayes  hyperprior 
approach  to  empirical  Bayes  problems  and  the  way  the  estimates  are 
influenced  by  the  prior  choice,  an  example  is  given  here. 

The  data  set  of  this  example  is  taken  from  Bayesian  Reliability 
Analysis  (p.  636)  by  Martz  and  Waller.  Suppose  a  10  hour  life  test 
has  been  conducted  for  each  of  the  eleven  production  lots  of  a  high 
reliability  device.  The  numbers  of  failures  were  observed  to  be 
0,  1 ,  0,  0,  1 ,  2 ,  0,  1 ,  0 ,  0,  and  0  respectively.  It  is  assumed 
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that  the  event  of  failures  can  be  modeled  by  Poisson  point  processes 

with  intensity  rate  for  the  ith  lot.  The  objective  is  to  estimate 

5 

for  each  lot.  Let.  6 ^  =  A ^  *  10  .  A  Dirichlet  hyperprior  approach 


to  this  problem  is  to  assume  are  distributed  according  to  an  unknown 

distribution  G  .  G  is  chosen  from  the  Dirichlet  process  with  prior 
&{(-°st]}  =  MGqU)  ,  where  Gq  represents  the  statistician's  prior 
guess  of  G  ,  and  M  represents  the  statistician's  strengths  of  prior 


belief  in  G^  .  Then  the  Bayes  estimator  of  can  be  obtained  from 


( 1 ) ,  wnere 


"ei  pxi 
e  0i 


fi(xi!ei}  =  —771 


x.  denotes  the  number  of  failures  in  the  10  hour  test  in  the  ith  lot. 

The  evaluation  of  the  Bayes  rule  (1)  can  be  approximated  by  an 
importance  sarrpling  Monte  Carlo  method  proposed  by  Kuo  (1985).  If  we 
choose  Gn  a  gamma  distribution  G(a,S),  i.e., 


G0(d; 


e“  cc.-i  -s- 

7 — r  6  e 


1(6  >  0)  , 


then  the  single  integrals  (see  equations  (£)  and  (5)  of  Kuo,  1985) 
contained  in  each  Monte  Carlo  iteration  can  be  evaluated  by  using 
the  followinG  identities: 


/  6  7,  f(xi  e-)GQ(dr)  = 
i  c  k 


ea  r(a  -  1  *  Tk  x,) 


T(a) • ( ! k  j  +  6) 


J1  *  1  +  iL  N> 


v  X  . 

ick  1 


/  r  f(xi le)G0(d5) 

ick 


e°  r(*  +  Tk  *,) 


r(o)-(!k|  *  e) 


h  *  uk  xi> 


.n  x 

ick  *i 
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where  k  is  a  subset  of  the  index  set  {1,  ....  11},  ]k|  denotes  the 

number  of  indices  in  k  . 

If  we  choose  Gq  a  uniform  distribution  on  (0,  0^),  then  the 
single  integrals  are  evaluated  by 


/  a  f(xJe)G0(dr)  =  llkV  2  *  Uk  V  r(;  *  kk  V  , 


e0  •  |k;'2  +  ilk  «,)  •  &  x. ! 


ick  i 


I  Uk  f(«,|5)s0(ae) 


=  Kks„,  l  *  A  xj  r(l  *  X  x) 


k|(1  +  ikxi> 


■  IT  I  x-  ! 
ick  i 


y  r_i 

where  I(y,  r)  =  /  t  e  d t/r ( r )  denotes  the  incomplete  gamma 
0 

function . 

In  the  following  tables,  is  chosen  to  be  either  a  gamma 

G(:i,£)  or  a  uniform  distribution  u(0,5q)  .  Tne  Bayes  rules  (1)  for 
various  values  of  a,  6,  0q  and  M  are  computed.  They  are 
evaluated  by  Monte  Carlo  methods  with  the  number  of  iterations 
NI  =  4000  or  NI  =  16000  .  Therefore,  the  posterior  standard  errors 
(see  equation  (12)  of  Kuo,  1985)  are  included  in  parentheses.  We 

have  also  rearranged  the  order  of  the  observations  for  easier  visual 

A  ~  5 

examination.  The  intensity  rate  A.,  are  estimated  by  \ ^  -  0^/10 

In  Table  1,  a,  6  for  the  gamma  distribution  prior  guess  are 
*  -2  2-  ~ 

chosen  to  be  a  =  max  {x  /(S  -  x),  0}  ,  and  8  =  max  {x/(S  -  x),  0} , 

where  x  =  [x^/11  ,  and  =  J(x.  -  x)  / 1 0  .  This  choice  of  prior 
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is  motivated  by  assuming  that  the  0.’s  are  independent  and  identically 
distributed  according  to  G( a, 8)  ,  where  a  and  B  are  unknown. 
Considering  the  marginal  distribution  of  x  ,  we  have  E  =  a/B  ,  and 

A  /S 

V(x.)  =  a(l  +  l/p)/3  "  for  all  i  .  The  estimators  a  and  8  above 
are  obtained  by  the  method  of  moments  and  adjusted  for  a  >  0  and 
6  >  0  . 

Robbins  (15SG)  discussed  two  methods  of  obtaining  empirical  Bayes 
estimators  and  proposed  a  method  for  combining  the  two  to  achieve  both 
consistency  and  efficiency.  The  first  estimator  for  6^  denoted  by 
6.  ,  is  derived  from  a  nonparametric  prior  point  of  view.  Then 
-•j(*)  =  (x.j  +  1)  *  #(x.  +  l)/#(x-)  ,  where  #(x)  denotes  the  number  of 
lots  with  x  items  failed.  The  second  estimator  for  0^  ,  denoted  by 
x .j  ,  is  derived  from  a  parametric  gamma  prior: 

t^x)  =  (Xj  +  a)/(  1  +  S)  =  x1  -  (x/S^ )  •  I(x  <  S«)  (x^  -  x)  , 

where  x.  >  0  for  all  i  .  These  are  Stein  type  shrinkage  estimators. 
Both  types  of  estimators  are  given  in  Table  1.  For  x^  ,  we  use  the 
unbiased  estimator  for  the  variance  of  the  x's  instead  of  the  usual 
sample  variance  estimator  suggested  by  Robbins.  Otherwise,  the 
shrinkages  of  the  data  will  pass  the  origin  without  the  positive  part 
corrections  given  in  the  x^  expression  above. 

For  the  Dirichlet  hyperprior  approach,  with  the  G^.  chosen  as 

/V  /V 

G(a,6)  above,  we  would  expect  the  proper  Bayes  estimators  to  approach 
to  the  x's  as  M  — >  °°  .  This  is  confirmed  by  the  estimators  shown 
in  Table  1 . 


10 


In  Tables  2(a)  and  2(b),  a  different  shape  of  the  gamr.a  distribution 

is  chosen.  The  mean  of  in  Table  2(a)  is  the  same  as  before,  i.e., 

x  .  However,  the  variance  of  Gq  is  4.54544  which  is  much  larger  than 

0.0181817,  the  one  used  in  Table  1.  It  is  interesting  to  observe  that 

the  estimators  are  closer  to  the  maximum  likelihood  estimator  as  M  — >  <*. 

This  reveals  that  choosing  the  prior  on  a  mixing  distribution  is  quite 

different  from  the  usual  nonparametric  prior  situation  without  the  mixing 

situation.  As  M  — >  °°  ,  we  expect  the  9’s  are  more  distinct.  When 

11 

M=5G0,  P(all  8's  are  distinct)  =  (1  -  ( i -1 )/(M+i -1 )) =  0.897. 

Therefore,  the  estimator  6^  can  be  approximated  by  using  x^  alone  and 
the  prior.  When  M  is  small,  many  of  the  0 ' s  are  identical  with  high 
probability.  This  explains  the  phenomenon  exhibited  in  Table  2(a)  and 
2(b). 

The  G^'s  in  Tables  3(a)  and  3(b)  have  the  same  mean  and  different 
variances.  Different  means  and  different  variances  are  selected  for 
Tables  4(a)  and  4(b)  to  examine  how  the  estimates  are  influenced  by  the 
prior  guess. 

In  Tables  5(a)  and  5(b),  the  Gq's  are  chosen  to  be  the  uniform 
distributions  U(0,0.3)  and  U(G,1)  respectively.  The  incomplete  gamr.a 
distributions  described  in  (2)  are  evaluated  by  the  MDGAM  subroutine  of 
the  International  Mathematical  and  Scientific  Library  (IMSL). 


n 


r  TABLE  1 

Bayes  and  Empirical  Bayes  Estimates 
for‘  {0.}  with  Go  ~  G(u,  £) 
a  =  11.3636  S  =  25  NI  =  4000 


X  . 

1 

^  (V.=  l) 

9.  (M=50) 

ei  (m*ioo) 

6 .  (M=500) 

C1 

Ti 

0 

0.4441  (0.0005) 

0.4372  (0. 0002) 

0.4372  (0.0001) 

0.4370  (0.0001) 

0.4286 

0.4371 

0 

0.4453  (0.0C05) 

0.4372  (0.0C02) 

0.4371  (0.0001) 

0.4370  (0.0001) 

0.4236 

0.4371 

0 

0.4449  (0.0005) 

0.4371  (0.0002) 

0.4371  (0.0001) 

0.4370  (0.0001) 

0.4286 

0.4371 

0 

0.4440  (0.0005) 

0.4373  (0.0002) 

0.4370  (0.0001) 

0.4369  (0.0001) 

0.4236 

0.4371 

0 

0.4453  (0.0005) 

0.4372  (0.0002) 

0.4371  (0.0001) 

0.4370  (0.0001) 

0.4286 

0.4371 

0 

0.4451  (0. 0005) 

0.4363  (0.0002) 

0.4369  (0.0001) 

0.4369  (0.0001) 

0.4286 

0.4371 

0 

0.4455  (0.QC05) 

0.4375  (0.0002) 

0.4372  (0.0001) 

0.4369  (0.0001) 

0.4286 

0.4371 

1 

0. 4635  (0.0005) 

0.4743  (0.0002) 

0.4749  (0.0001) 

0.4754  (0.0001) 

0.6667 

0.4755 

1 

0.4644  (0.0005) 

0.4749  (0.0002) 

0.4753  (0.0001) 

0.4753  (0. 0001 ) 

0 . 6667 

0.4755 

1 

0.4549  (0.GC05) 

0. 47d2  (0.0002) 

0.4751  (0.0001) 

0.4754  (0.0001) 

0.6667 

0.4755 

2 

0.4356  (0.0004) 

0.5125  (Q.G001) 

0.5132  (0.0001) 

0.5133  (0.0000) 

0.0000 

0.5140 

12 


TA3LE  2 


Bayes  Estimates  for  {?.} 
"  with  ~  G(a,  S) 

(a)  a  =  0 . : A- 4544  3  =  0.1  M  =  15000 


X  . 

1 

2  1 

;v=d 

5*  (M=5G) 

6  (M= 100) 

0  (M=500) 

0 

0.2:30 

(C.01C;) 

0.0:90  (0.0C17) 

0.C512  (0.0012) 

0.0430  (0.CCC2) 

0 

0.22-32 

(u. 0097) 

0.0597  (0.0C19) 

0.0510  (0.0007) 

0.C430  (0.0002) 

0 

C.2925 

(0. Cl  02) 

O.Cc1!  (C. 0021 ) 

0.0499  (0.0006) 

0.0433  (0.0002) 

0 

0 . 2  '  0  5 

(0.  Ct'33) 

0.C6C7  (0.CC19) 

0.0493  (0.0CC6) 

0.0433  (0.G002) 

0 

0.2233 

(0.0101) 

0.0636  (0.0017) 

0.0513  (0.0012) 

0.0430  (0.0002) 

a 

O.icCl 

(0.033c) 

0.0576  (0.0015) 

0.0511  (0. 0012) 

0.0430  (0.C002) 

0 

0.2:55 

(0.0102) 

0.0353  (0.0014) 

0.0514  (0.0007) 

0.0429  (0. 0002 ) 

l 

0.7376 

(0.0123) 

0.9771  (0.0067) 

0.9624  (0.0024) 

0.9573  (0.0019) 

l 

0. 7554 

(G. 0129) 

0.9715  (0.CC69) 

0.9647  (0.0039) 

0.9534  (0.0011) 

l 

C. 7307 

(0.0130) 

0.9751  (0.0067) 

0.9617  (0.0023) 

0.9545  (0.0012) 

2 

0. 31 54 

(C. Clio) 

1.6952  (0.0111) 

1.7753  (0.0040) 

1.8329  (0.0020) 

(b) 

a  =  0.1 

3  =  0.01 

NI  =  16000 

x . 

l 

5.  (»1 ) 

6  ( .V=  50) 

0  (M=100) 

0  (M=530) 

0 

0.3131  (0.CC31) 

0.1155  (0.G012) 

0.1076  (0.0003) 

0.1004  (0.0002) 

0 

0.3227  (0.0031) 

0.1163  (0.0013) 

0.1084  (0.0C07) 

0.1004  (0.0003) 

0 

0.3243  (0.0032) 

0.1175  (0.0014) 

0.1070  (0. 0006) 

0.1006  (0.0003) 

0 

0.3133  (0.0081) 

0.1170  (0.0012) 

0.1070  (0.0006) 

0.1007  (0.0003) 

0 

0.3253  (0.0034) 

0.1174  (0.0012) 

0.1030  (0.0008) 

0.1003  (0.C003) 

0 

0.3194  (0.0079) 

0.1'43  (0.0011) 

0.1077  (0.0003) 

0.1004  (0.0003) 

0 

0.3263  (0.0081) 

0.1138  (0.0011) 

0. 1033  (0.0007) 

0.1002  (0.0002) 

1 

0.7163  (0.0115) 

1.0595  (0.0026) 

1.0733  (0.0016) 

1.0360  (0.0009) 

1 

0.7123  (0.0116) 

1.056J  (0.0027) 

1.0743  (0.0017) 

1 .0344  (0.0007) 

1 

0.7073  (0.0116) 

1.0535  (0.0026) 

1.0715  (0.0017) 

1.0550  (0.0008) 

i 

0.8632  (0.0094) 

1.9243  (O.uCSl ) 

1.9947  (0.0029) 

2.0552  (0.0014) 

n 


TABLE  3 


Dd_.es  Estimates  fcr  {5} 
with  G0  "  G{i,  £) 


(a) 

-  : 

-7 - 

3  =  50 

Nl 

:  =  4, 

330 

Xi 

'  i 

h  < 

v  i 

(7- ICC) 

rj 

V 

G. 

--.7  3 

/  -V  ^  \ 

0. 

fl  .1  ^  £ 

.  :o 

(O.CQCl ) 

0. 

443  5 

(0.0001) 

W 

0. 

1  ~  -  “ 

(3.0:03; 

0. 

.  **  -f  5  / 

(0.C001) 

0. 

442  5 

(0.0001) 

0 

G. 

4:3  > 

(0.0:03 ) 

0, 

.445  7 

(0 . 0001 ) 

0. 

4455 

(O.CCOl) 

c 

0. 

4  4  37 

(0 . c:c3) 

0. 

,4455 

(0.0301) 

0. 

4455 

(O.COOl) 

G 

0. 

(o.c:c3) 

0. 

.4455 

(O.COOl) 

0. 

44  35 

(C.C331) 

0 

0. 

4^9: 

(C.COOO) 

0. 

*  ■ 

,  *t  t  _  J 

(0.C031) 

0. 

4455 

(O.CCOl) 

0 

0, 

4oGG 

(0.0003) 

0. 

,4455 

(O.CCOl) 

0. 

4455 

(O.COOl) 

1 

0. 

4593 

(0.0003) 

0. 

.4549 

(0.0001) 

0. 

4649 

(G.0001) 

1 

0. 

4593 

(0.C003) 

0. 

.4647 

(0.0001) 

0. 

4651 

(O.COOl) 

1 

0. 

4501 

(0.G0C3) 

0. 

.4549 

(O.COOl) 

0. 

4650 

(O.CCOl) 

2 

G. 

4c33 

(0.0332) 

0. 

,4241 

(O.CCOl) 

0. 

4344 

(0.00C1) 

(A) 

}  : 

=  5.51 

^  1  : 

D  =  1 

2.5 

m 

.  =  40 CO 

xi 

9  ; 

1 

;m*i) 

( 

7=50) 

9i 

(M=i.:o) 

0 

0. 

.4354 

(d.c:c?; 

0, 

.4213 

(0.0003) 

0. 

4213 

(O.OCC2) 

r* 

U 

0. 

4372 

(0  .CO/:*) 

0. 

,4215 

(0.GCC3) 

0. 

4212 

(O.COC2) 

0 

0 . 

.  437i 

(0.0009) 

0, 

,4215 

(0.0C03) 

0. 

4211 

(0.0002) 

0 

0. 

4352 

(o.cc 

0. 

4210 

(0.0003) 

0. 

4210 

(0.0002) 

G 

0. 

.4353 

(0.0009) 

0, 

.4210 

(0.0003) 

0. 

4210 

(0.CG02) 

G 

0, 

.4253 

(0.CCC9) 

0. 

,4211 

(0.0003) 

0. 

4203 

(0. 0002) 

0 

0, 

4373 

(0.0009) 

0, 

,4211 

(0.0C03) 

0. 

4213 

(0.CC02) 

1 

G, 

,4746 

(0.0010) 

0. 

4941 

(0.0003) 

0. 

4940 

(0.0002) 

1 

0, 

.4744 

( 0.CCC9) 

0, 

4935 

(G.0C03) 

0. 

49 '7 

(0.0003) 

1 

0, 

.4757 

(0.0010) 

0, 

,4941 

(O.CC03) 

0. 

4943 

(0. 0002 ) 

2 

G. 

,5091 

(0.0009) 

0, 

,  56c-2 

(0.0002) 

0. 

5674 

(0.0002) 
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TABLE  4 


Bayes  Estimates  for  {9-} 
with  Go  ""  G(cx,  B) 


(a) 

a  =  1 1  .3636 

8  =  12.5 

NI 

=  4000 

xi 

A 

6i  ( 

;m=i) 

A 

6i  ' 

;m=5o) 

ei  (m=ioo) 

0 

0.7493 

(0.0012) 

0.3366 

(0.0004) 

0.8394 

(0.0003) 

0 

0.7502 

(0.0012) 

0.8371 

(0.0004) 

0.8392 

(0.0003) 

0 

0.7500 

(0.0012) 

0.8366 

(0.0004) 

0.8391 

(0.0003) 

0 

0. 7486 

(0.0012) 

0.8358 

(0.0004) 

0.8389 

(0.0003) 

0 

0.7496 

(0.0012) 

0.8361 

(0.0004) 

0.8391 

(0.0003) 

0 

0.7495 

(0.0017) 

0.8364 

(0.0004) 

0.8389 

(0.0003) 

0 

0.7501 

(0.0012) 

0.8362 

(0.0004) 

0.8391 

(0.0003) 

1 

0.7851 

(0.0015) 

0.9093 

(0.0004) 

0.9120 

(0.0003) 

1 

0.7846 

(0.0015) 

0.9085 

(0.0004) 

0.9128 

(0.0003) 

1 

0.7871 

(0.0015) 

0.9091 

(0.0004) 

0.9122 

(0.0003) 

2 

0.8232 

(0.0018) 

0.9816 

(0.0004) 

0.9857 

(0.0003) 

(b) 

X  ■ 
l 

a  =  10 

A 

Cl  l 

"  1 

(M=l) 

B  =  50 

A 

9i  ' 

NI 

;m=5o) 

=  4000 

6i  ( 

;m=ioo) 

0 

0.2224 

(0.0004) 

0.1972 

(0.0001) 

0. 1968 

(0.0001) 

0 

0.2227 

(0.0004) 

0.1973 

(0.0001) 

0.1968 

(0.0001) 

0 

0.2229 

(0.0004) 

0.1974 

(0.0001) 

0.1968 

(0.0001) 

0 

0.2225 

(u . 0004) 

0.1972 

(0.0001) 

0.1967 

(0.0001) 

0 

0.2226 

(0.0004) 

0.1972 

(0.0001) 

0.1967 

(0.0001) 

G 

0.2221 

(0.0004) 

0.1972 

(0.0001) 

0.1966 

(0.0001) 

0 

0.2227 

(0.0004) 

0.1972 

(0. 0001 ) 

0.1968 

(0.0001) 

1 

0.2331 

(0.0003) 

C. 21 66 

(0.0001) 

0.2161 

(0.0001) 

1 

0.2333 

(0.0003) 

0.2165 

(0.0001) 

0.2163 

(0.0001) 

1 

0.2333 

(0.0003) 

0.2167 

(0.0001) 

0.2162 

(0.0001) 

2 

0.2422 

(0.0002) 

0.2358 

(0.0001) 

0.2356 

(0.0001) 
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TABLE  5 

Bayes  Estimates  for  {9^} 
with  G0  ~  U(0t  e0) 


(a)  90  =  0.3  ,  NI  =  4000 


xi 

6.  (M=1 ) 

0.  (M=50) 

0.  (M=100) 

0 

0.2003  (0.0003) 

0.1455  (0.0002) 

0.1442  (0.0002) 

0 

0.2007  (0.0008) 

0.1457  (0.0002) 

0.1442  (0.0002) 

0 

0.2013  (0.0007) 

0.1458  (0.0002) 

0.1442  (0.0002) 

0 

0.2007  (0.0007) 

0.1456  (0.0002) 

0.1441  (0.0002) 

0 

0.2005  (0.0008) 

0.1455  (0.0002) 

0.1441  (0.0002) 

0 

0.1997  (0.0008) 

0.1454  (0.0002) 

0.1439  (0.0002) 

0 

0.2008  (0.0007) 

0.1455  (0.0002) 

0.1443  (0.0002) 

1 

0.2222  (0.0004) 

0.1965  (0.0002) 

0.1956  (0.0001) 

1 

0.2225  (0.0005) 

0.1963  (0.0002) 

0.1958  (0.0001) 

1 

0.2224  (0.0004) 

0.1966  (0.0002) 

0.1957  (0.0001) 

2 

0.2323  (0.0002) 

0.2223  (0.0001) 

0.2219  (0.0001) 

(b) 

e o  *  i.o  hi  » 

16000 

x . 

1 

ei  (M=1) 

0i  (M=50) 

0.  (M=7 00) 

0 

0.4601  (0.0015) 

0.4192  (0.0004) 

0.4185  (0.0003) 

0 

0.4627  (0.0014) 

0.4195  (0.0004) 

0.4183  (0.0003) 

0 

0.4635  (0.0015) 

0.4194  (0.0004) 

0.4183  (0.0003) 

0 

0.4615  (0.0015) 

0.4195  (0.0004) 

0.4185  (0.0003) 

0 

0.4635  (0.0014) 

0.4194  (0.0004) 

0.4182  (0.0003) 

0 

0.4625  (0.0014) 

0.4189  (0.0004) 

0.4183  (0.0003) 

0 

0.4635  (0.0014) 

0.4185  (0.0004) 

0.4187  (0.0003) 

1 

0.5679  (0.0011) 

0.6059  (0.0003) 

0.6067  (0.0002) 

1 

0.5673  (0.0011) 

0.6056  (0.0003) 

0.6069  (0.0002) 

1 

0.5661  (0.0011) 

0.6059  (0.0003) 

0.6066  (0.0002) 

2 

0.6301  (0.0008) 

0.7058  (0.0001) 

0.7072  (0.0001) 
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